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Abstract 

A previous analysis of the charge (Z) correlations in the AZ— < Z > plane for 
^^^Xe+""*Sn central collisions at 32 MeV/n has shown an enhancement in the pro- 
duction of equally sized fragments (low AZ) which was interpreted as an evidence 
for spinodal decomposition. However the signal is weak and rises the question of the 
estimation of the uncorrelated yield. After a critical analysis of its robustness, we 
propose in this paper a new technique to build the uncorrelated yield in the charge 
correlation function. The application of this method to ^^^Xe-|-""*Sn central colli- 
sion data at 32, 39, 45 and 50 MeV/n does not show any particular enhancement 
of the correlation function in any AZ bin. 
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1 Introduction 



The search for signals of a liquid-gas phase transition in the nuclear medium 
is presently a very active field of research (see [1] and ref. therein). In such 
a framework, multifragmentation is a key process. A possible scenario is the 
excursion of the excited system, produced in dissipative heavy-ion collisions, 
at low densities inside the coexistence region of the nuclear phase diagram (for 
a general introduction to the subject, see for instance [2]). There, a volume 
instability, the so-called spinodal decomposition, leads to a disassembly of the 
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system [3,4]. Such a behaviour is predicted by microscopic transport calcula- 
tions based on the nuclear Boltzmann equation and its stochastic extensions 
[5,6,7]. The main idea is to search for fragmentation events in which most 
fragment atomic numbers are nearly equal as suggested by theoretical works 
[8,9]. It is worth noting that surface instabilities, like Rayleigh or sheet insta- 
bilities [10,11], may produce, in the early compression-expansion stage of the 
collision, the formation of toroidal and bubble nuclei (depending of the value 
of the incompressibility constant K) that may also decay into several nearly 
equal fragments [12]. In any case, the proportion of this kind of multifragmen- 
tation events should be small and a sensitive method, called "higher-order 
correlations", has been proposed in Ref. [13] to search for weak signals in 
the fragment charge correlation analyses. Following the prescriptions of this 
method, the question of a possible experimental signature of spinodal decom- 
position of excited nuclear matter has been examined for ^^^Xe-|-""*Sn central 
collisions at 32 MeV/n in a previous work [14]. 

The analysis relies on two variables. The average fragment charge per event is 
defined as: 

1 M 

<z>=—yzi (1) 



where Zj is the charge of fragment i and M is the fragment multiplicity. Only 
fragments with atomic number Zj > 5 are considered. As in Refs [13,14], the 
standard deviation in fragment sizes is defined as: 
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The correlation function, i?, is built by considering the ratio of the measured 
yield F(AZ, < Z >) versus the "uncorrelated" yield F' (AZ, < Z >). Since 
events with nearly equal size fragments are associated with low values of AZ, 
any significant enhancement of R for small AZ (AZ < 1) may signal a spin- 
odal decomposition as proposed in [14]. The key question is how to build the 
uncorrelated yield Y' . This last quantity has been estimated for each class of 
fragment multiplicities by mixing fragments from different events of the se- 
lected sample [13,14]. In doing so, the initial total charge carried out by the 
fragments {Zhound = Y^fii Zi = M < Z >) and also the average fragment 
charge < Z > are not conserved in building Y'. Consequently, in the absence 
of correlations the ratio R = Y/Y' is not equal to unity and depends simul- 
taneously on the < Z > and AZ values. Thus, to extract the true correlation 
one has to determine a background function R' corresponding to the ratio 
Y/Y' for uncorrelated events. 
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In the present paper, in a first step, we would like to show how the method 
previously used to generate the charge correlation functions [14], that we will 
call Zhound mixing method, may generate spurious signals. Then, a new method 
to build the uncorrected yield Y' is presented in Section 3. By applying this 
method, the fragment charge correlations obtained in the ^^^Xe+"°*Sn central 
collisions between 32 and 50 MeV/n are discussed in Section 4. A recently 
proposed method to build the denominator Y', called IPM for Intrinsic Prob- 
ability Method [15], which has been applied to the same Xe+Sn data set [16], 
is commented in Section 5. 



2 Review of the Zhound mixing method 

2. 1 Data selection and previous charge correlation results 

The charge correlations were analysed in a set of events associated with a 
"quasi-fusion" process occurring in central ^^sXe+^'^^Sn collisions at 32 MeV/n. 
Only events with the sum of the detected charges greater or equal to 80 
were considered, the total charge of the Xe+Sn system being 104. Briefly, 
the event selection was based on the measurement of the distribution of the 
preferred emission direction of matter in the centre of mass of the reaction 
(flow angle). In this work, fragment charge correlation analyses have been 
performed by selecting events having a flow angle larger than 60°. It has been 
shown [17,18] that such events {Qfiow > 60°) could be interpreted as resulting 
from the multifragmentation of a thermalised and compact nuclear system 
("single source" events). In a previous work [14], using the same gating condi- 
tion, the charge correlation functions R{AZ, < Z >) were first built for each 
fragment multiplicity M = 3, 4, 5 and 6. However, the results suffered from a 
lack of statistics. The numbers of events in the first bin in AZ {AZ < 1) and 
for each multiplicity of this selection are shown in the third column of Table 1. 
As we see, these numbers, summed over all < Z >, are quite small, specially 
for multiplicities M = 5 and 6. Thus, to increase the statistical significance of 
the signal, the charge correlations were built by replacing the variable < Z > 
by the sum of the fragment charges Zhound (each bin in Zhound containing six 
atomic numbers). Summing over all M, the resulting correlation functions 
shown in Fig. 3 of Ref. [14] were interpreted as an evidence of a significant 
enhancement of equal-sized fragments associated with fused events in Xe+Sn 
collisions at 32 MeV/n. However, as we mentioned previously, the charge cor- 
relation function cannot be compared to unity but to a specific R' background 
function. This function has been estimated in Ref. [14] by assuming that, in 
the absence of correlations and for a given Zbound, R' = f{AZ, Z^ound) is expo- 
nentially increasing with AZ in the range < AZ < 8. Obviously, the choice 
of the background estimation is the critical point of this method and in or- 
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der to test the hypothesis of Ref. [14] we have performed a charge correlation 
analysis R{AZ, Zi,ound) from a sample of randomly generated events. 



2.2 Random break-up simulation 

In the Monte Carlo simulation, the following procedure, that we will call ran- 
dom break-up, has been used: for a given fragment multiplicity M, the charge 
of the fragments is obtained by breaking M — 1 bonds at random in a chain of 
length Zhound- Only events for which all fragments have Z > 5 were kept. For 
each multiplicity, the Zhound distribution was the experimental one for central 
Xe-|-Sn collisions at 32 MeV/n. We have checked that for each multiplicity M, 
the ensemble thus generated is equivalent to the one obtained by computing 
all partitions and their permutations containing M fragments with Zi > 5 
such that J2iiiZi = Zhound [19]. The Zhound distributions are given in the 
first column of Fig. 1 as a function of the fragment multiplicity M . They are 
gaussian in shape with a centroid which increases from ~54 for M=3 to ^64 
for M=6, whereas their width is practically constant. To minimize statistical 
effects, twenty times more events were generated in the simulation than in the 
data. The simulated Z-distributions (full lines) are compared to the experi- 
mental data in the second column of Fig. 1. For multiplicities M= 5 and 6, 
the agreement is remarkably good. At lower multiplicities M= 3 and 4, the 
simulated distributions are slightly broader than the experimental ones. The 
last two columns of Fig. 1 compare the simulated AZ distributions (full lines) 
to the experimental ones with two AZ scales (in order to show the details 
at low AZ). There again for multiplicities 5 and 6, the agreement between 
the simulation and the data is remarkable. For the lower multiplicities 3 and 
4, as for the Z-distributions, the simulated distributions are broader than 
the experimental ones. Examination of the last column of Fig. 1 reveals no 
enhancement over the simulation for AZ < 6. The numbers of events with 
AZ < 1 obtained in the simulation are listed in the second column of Table 1 
as a function of multiplicity. We note that the total fraction of events (0.154 
± 0.005 %) in this AZ bin is compatible with the experimental one (0.127 ± 
0.020 %). 

Fig. 2a shows the ratio between the random break-up yield Ybreak-up over 
ypartitions obtained by calculating analytically all possible partitions [19] as 
a function of AZ for each multiplicity M and their sum. Except for the low 
statistics data with large AZ (M=6), all of the ratio values are consistent with 
unity. It is worth noting that in the present case, deviations from unity are of 
statistical origin only. Fig.2b displays the values of Ybreak-up/ ypartitions in the 
bin AZ < 1 as a fonction of < Z >. In this bin, the statistics is poor and we 
observe deviations in excess of one standard deviation (a). In the following, 
the random break-up events are analysed using the Zhound mixing method [14]. 
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2.3 Application of the Zbound mixing method to the random break-up simula- 
tion 

Firstly, the charge correlation functions R{AZ, < Z >) were built for each 
multiplicity M = 3,4, 5, 6, separately. For this purpose, the uncorrelated yield 
Y' was built generating 10'^ times more events than in the numerator in order 
to minimize statistical errors on the denominator. The R values corresponding 
to four (M, < Z >) pairs chosen in the interval 48 < Zfyound < 52 are shown in 
Fig. 3a. They illustrate the general evolution of the function R = Y/Y'. The 
amplitude and extension in AZ of the function R depend on the multiplicity 
M and the average charge < Z >. However, we note that within statistics, 
the charge correlation functions decrease with decreasing AZ without any 
significant enhancement in the first bin in AZ, as expected for uncorrelated 
events. In the absence of correlations other than charge conservation, these 
results (Fig. 3a) represent the background function R' discussed above. Large 
fluctuations and statistical error bars are present in the low AZ bins. Thus, a 
correct estimate of the function R' from experimental data, often associated 
with weaker statistics, appears to be a very difficult task. 

Secondly, following [14], we built the charge correlations R{AZ, Zbound)- They 
are displayed in Fig. 3b for six values of Zbound = 39,45,51,57,63,69 as a 
function of AZ (each bin in Zbound regrouping six atomic numbers). The four 
(M, < Z >) pairs shown in Fig. 3a are now integrated in the bin Zbound = 51. 
Obviously, our events being uncorrelated all these correlation functions are the 
background functions R' . Contrary to what was assumed in [14] these back- 
ground functions are not exponential functions of AZ even in the restricted 
range AZ < 8. Thus, if as in [14], one assumes an exponential background in 
the range 2 < AZ < 8 (the lines in Fig. 3b), and extrapolates it down to zero, 
one creates artificially a signal at low AZ. The fragment charge correlations 
for the bin AZ < 1 and the corresponding background extracted this way 
(broken line) are plotted in Fig. 4 (up-left) as a function of Zbound- This figure 
is very similar to Fig. 3 of Ref. [14]. This spurious result finds its origin in 
the mixing of events with different fragment multiplicities M = 3, 4, 5 and 
6 to build Zbound bins. Indeed, the charge correlation functions R{AZ, Zbound) 
combine various functions R{AZ, < Z >) which are multiplicity dependent. 

This artificial enhancement of equal-sized fragments is not systematic and 
depends on the Zbound distributions. Xe-|-Sn central collisions have also been 
studied with INDRA at 39, 45 and 50 MeV/n [20,21]. In Ref. [21], the au- 
thors show the disappearence of the charge correlation signal at low AZ when 
the energy increases. But, the same evolution is obtained by analysing simu- 
lated data without charge correlations. By using the Zbound mixing method, 
we show in Fig. 4 the charge correlations in the first AZ bin of four sets of 
events generated by the random break-up simulations obtained from the Zbound 
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distributions of Xe+Sn data at 32, 39, 45 and 50 MeV/n (with 20 times the 
statistics of the data). This figure, very similar to Fig.2 of Ref. [21], shows 
clearly the artificial decreasing of the charge correlation signal when the en- 
ergy increases. Consequently, the method using the mixing of multiplicities 
[14] renders unreliable an estimation of the background R' . 



3 The exchange mixing event method 

3. 1 Description of the method 

To get rid of the problematic evaluation of the background, we attempted to 
build an uncorrelated yield Y' with the conservation of the total fragment 
charge [Zbound)- In such a case the charge correlation function can be directly 
compared to unity. The following procedure has been developed. For an ensem- 
ble of events with a given fragment multiplicity, M, we consider the exchange 
of two fragments of charges Zi{i) and Z2{i) belonging to the event labelled i 
with two fragments of charges Z[{j) and Zl^lJ) belonging to the event labelled 
j under the condition that Zi{i) + Z2{i) = Z[{j) + Zi^i^j). As such, the total 
fragment charge Zbound and the average charge < Z > are conserved event 
by event and then the Zbound and the charge distributions of the considered 
ensemble are conserved. For a given event we performed the substitutions 
by a random picking of event j followed by a random picking of charges inside, 
respectively, the events i and j. If the condition Zi{i) + Z2{i) = Z[{j) + Zi^^j) 
was verified the charge substitutions were performed, if not we proceeded with 
a new random picking of event j, and so on, until all events in the sample has 
been changed at least once (one iteration). First, 100 iterations are performed 
to get a mixed sample. Then, the yield F' (AZ, < Z >) is created from this 
mixed sample by continuing the same iterative procedure 10'^ times and at 
each time Y' is incremented in order to make negligible its statistical error 
(relative to Y). Accurate tests showing the sensitivity of the exchange method 
are developed in the next. One can also consider exchanges among the events 
in which more than two (say N) fragments are exchanged. In this latter case, 
the method is obviously restricted to events for which M is larger than and 
in the following we only use N = 2. 

3.2 Tests of the exchange mixing method 

The technique has been tested with the random break-up calculation and the 
results are plotted in Fig. 5a for each multiplicity and the whole ensemble 
M = 3,4,5,6. For any multiplicity, the correlation function R = f[AZ) is 
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constant and close to unity (most of the statistical error bars are smaller than 
the symbol size). Although we performed this simulation with a statistics 
twenty times larger than the one of the Xe+Sn data at 32 MeV/n (see Table 
1), statistical fluctuations are still present, specially, in the first and largest 
bins in AZ for M = 6. The correlation function R corresponding to the whole 
multiplicity set is equal to unity (bottom picture in Fig. 5a) independently 
of AZ. In this latter case, the fact that the correlations R = f{AZ) are 
multiplicity independent validates the multiplicity regrouping. 

However, when this method is applied to the data, the denominator may par- 
tially keep the fragment charge correlations present in the experimental sample 
(numerator). For example, applied on an event sample populated in majority 
by nearly-equal sized fragments, this method, as any others based on the data, 
would paradoxally give a weak signal. So, the exchange method may only be 
relevant in the search of a weak proportion of equal-sized fragment events in 
the data and it is crucial to test its sensitivity to measure an enhancement 
of charge correlations precisely in the bin AZ < 1. In order to do this test, 
we doubled the number of events populating the first bin in AZ in the ran- 
dom break-up sample and build new uncorrelated yields Y'. The results of the 
charge correlation analyses of this new sample are displayed in Fig. 5b. As 
expected, the charge correlations R in the first bin, within statistical uncer- 
tainties, are twice the ones shown in Fig. 5a, whereas they are unchanged in 
other bins. By integrating over all multiplicities M and < Z >, the correlation 
value in the first bin is: R = 1.941 ±0.045 (bottom picture in Fig. 5b). In this 
case, the number of events is large and we performed the same analysis with 
a set of simulated events statistically compatible with the Xe+Sn data at 32 
MeV/u. In the first bin, the new correlation value R = 1.961 ± 0.192 is very 
close to the expected value = 2 at ~ 5 standard deviations from one. These 
tests substantiate the good sensitivity of the exchange procedure used to build 
Y'. 

The use of the random break-up simulation is also a powerful tool to evaluate 
the kind of errors induced by the exchange mixing method. In such a simula- 
tion, where no correlation is expected, the deviations from unity should contain 
a statistical part due to the number of generated events (numerator), called 
cT/vT, and, perhaps, a systematic error due to the method itself. In Fig. 6a, we 
show the deviations from 1 of all charge correlation values R{AZ, < Z >, M), 
normalized to their statistical error bar aN (absissa x = {R — l)/aN), obtained 
from the simulations previously described and for each multiplicity M = 3, 4, 5 
and 6 (see Fig. 5a). The good superimposition of the x distribution with the 
normal function (see Fig. 6a) shows clearly that the errors associated with the 
exchange method are purely statistical. As all possible partitions (YpartiUons) 
have been analytically calculated (see Fig. 2) in the random break-up sim- 
ulation, the yield Y' = YpartiUons has been used as denominator in Fig. 6b. 
The remarkable similarity between both x distributions in Fig. 6a and 6b in- 
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dicates that there is no contribution due to the exchange method itself in the 
fluctuations of R. In the same manner, we plotted in Fig. 6c, the distribution 
of X obtained from the Xe+Sn data at 32 MeV/n. In this data, the good su- 
perimposition of the distribution of R{AZ, < Z >,M) values by the normal 
function, indicates that all deviations of R from unity appear compatible with 
statistical fluctuations. In the next Section, the fragment charge correlation 
functions associated to the Xe+Sn data are presented as a function of AZ, 
< Z > and M. 



4 Application of the exchange mixing method to the Xe+Sn data 

In the Xe+Sn central collisions, the multifragmentation process is present 
at 32, 39, 45 and 50 MeV/n and complete events {Ztot > 80) selected with 
Qfiow > 60° have been considered to arise from the multifragmentation of 
single nuclear sources [17,20]. For a given fragment multiplicity, the number 
of events with AZ < 1 at respectively 39, 45 and 50 MeV/n bombarding 
energies are listed in columns 4 to 6 of Table 1. The production of near equal- 
sized fragment events increases with the incident energy and reaches ~ 2% at 
50 MeV/n. This trend is a result of charge conservation. The multiplicity of 
light charged particles and fragments {Z < 4) is increasing with the incident 
energy [22]. Hence, the average total fragment charge (Zhound) decreases. Thus 
for a given fragment multiplicity, the size of the fragments decreases and the 
probability to measure equal-sized fragment events naturally increases. The 
analyses of the charge correlations for these "single source" samples have been 
performed, using the exchange mixing procedure to generate Y'. The charge 
correlations R have been built at each energy. First, we show in Fig. 7 the 
charge correlation functions R = Y/Y', summed over all < Z >, as a function 
of AZ for each multiplicity M = 3, 4, 5 and 6 (the first four rows) and for 
each incident energy (columns). The bottom row of Fig. 7 corresponds, at 
each bombarding energy, to the whole experimental data (summed over M) 
and benefits from enhanced statistics. All values of R are close to unity. To 
be more quantitative. Fig. 8 displays the quantity (i? — l)/cr as a function of 
AZ for random break-up (with analytical background) simulated events (first 
column) and for the experimental data from 32 to 50 MeV/n (columns 2 to 5), 
for multiplicities M = 3 to 6 (rows 1 to 4) and for all multiplicities (last row). 
For convenience, we have drawn dotted lines corresponding, respectively, to 
one and two standard deviations. We note that the fluctuations around zero 
are about of the same amplitude in the simulation and in the data over the 
whole range in AZ. The fluctuations for small AZ (as well as for large AZ) 
are not particularly larger than in any other AZ bins. Nevertheless, following 
theoretical predictions [8,9] the overproduction of equal-sized fragments should 
specifically concern a narrow range in Z ~ 10 — 15. If true, the summing over 
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all < Z > to determine the correlation functions may reduce significantly the 
importance of the "signal". For this reason, we show in Fig. 9 the R values 
corresponding to this first bin as a function oi < Z > for each multiplicity 
M = 3, 4, 5 and 6 (the first four rows), for all multiplicities (last row) and 
for each incident energy (columns). All < Z > bins populated with more than 
one event have been plotted in Fig. 9. This choice induces the absence of data 
point at 32 MeV/n for M =5 and 6. Notice that experimental error bars are 
statistical only. All values of R are within 2a of 1. It has to be emphasized 
that the R values (Fig. 9) for each multiplicity corresponding to the energies 
45 and 50 MeV/n are particularly close to the R = 1 line and associated 
with small statistical fluctuations. On the other hand, at 32 MeV/n and 39 
MeV/n, we are dealing with large statistical errors. The, perhaps, systematic 
enhancements of R at 39 MeV/n for multiplicity 4 in the range < Z >= 7 — 11, 
are compatible with 1 within 1.2a which is poorly significant (less than 80% 
of confidence level). Finally, we do not observe any significant enhancement 
of equal-sized fragments in the Xe+Sn central collisions between 32 and 50 
MeV/n. 



5 Comparison between methods 

Recently, an alternative method, called IPM [15], has been proposed to build 
the denominator Y' using the independent emission hypothesis constrained 
by the charge conservation. Basically, it is assumed that the probability to 
observe a given partition (n^ : (ni, rizt^t)) with a multiplicity M of a nucleus 
of charge Ztot is given by the multinomial formula 

(P )'^^ 

P({n4) = aM\ n with E = 1 (3) 

where a is the normalization constant {J2n^ f{{''^z}) = 1), Pz is the intrinsic 
probability for emitting the charge z, the number of charges z in the parti- 
tion {Ztot = Hz "^zz) and M the multiplicity of the partition (M = Y.z ^z)- We 
note that in the case of a fully random break-up, the partitions are weighted 
by the factor M!/ H'^^!- Thus, in that particular case the quantity rK-Pz)""" is 
a constant. For large Ztot, the normalization condition implies Pz = 2~^ [15]. 
Any other laws for the P^ will imply a departure from random break-up and 
the partitions thus generated might not be free from correlations. In the IPM 
method, the probabilities P^ are determined by fitting the partition probabil- 
ities P{{nz}) of the numerator and, as a consequence, they may keep memory 
of the experimental correlations. In other words, the IPM method may mea- 
sure charge correlation deviations from the independent emission hypothesis 
but the absence of charge correlations in the denominator Y' is not obvi- 
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ous. Moreover, as the IPM method rehes upon a fit to the statistical weights 
associated with each fragment partition, it washes partly out statistical fluc- 
tuations. In Ref. [16], the IPM method has been applied to the Xe+Sn data 
and the deviations of R from unity are distributed according to a law which 
is narrower than a normal distribution. However what should be this law is 
unknown. Thus, the fact that the fluctuations of R around unity are narrower 
in the case of the IPM method compared with the exchange mixing method 
does not guarantee that the former method is more sensitive. The question of 
sensitivity has been accurately discussed for the exchange mixing method by 
performing several tests from simulated event samples in which the amplitudes 
of charge correlations in the bin AZ < 1 were precisely known (see Section 
3.2). 

In [15,16], the IPM method has been applied to a set of partitions generated by 
the Brownian One-Body (BOB) dynamical model which supports the spinodal 
decomposition of hot and dilute finite nuclear systems [5,7]. This calculation, 
performed for head-on Xe-|-Sn collisions at 32 MeV/n, has been shown to 
reproduce quite well the experimental fragment multiplicity and charge dis- 
tributions [23,24]. Of the 1% of events in the bin AZ < 1, a percentage of 
0.36% of events has been evaluated in excess of the background with the IPM 
method which represents a correlation value i? ~ 1.5 [16]. The charge corre- 
lation function of the same BOB event sample has been also determined by 
using the exchange mixing method proposed in this work and, inside the sta- 
tistical error bars, no enhancement of equal-sized fragments was visible [25]. 
This contradictory result is really puzzling. Indeed, from the simulated event 
sample obtained from the Zi^ound distribution of the Xe-fSn data at 32 MeV/n 
(see Section 3.2), with a statistics comparable to the BOB calculation one, we 
increased the number of events in the first AZ bin by 50%. We analysed, by 
taking the definition AZ = I]j(-Zj— < Z >Y used in Ref. [16], the frag- 
ment charge correlations of this new set of events with the exchange method. 
In the first AZ bin the charge correlation value, integrated over all M and 
< Z >, was: R = 1.450 ± 0.121, which is compatible with the expected value 
R = 1.5 at 3.7a from R = 1. Once more, we conclude about the good abil- 
ity of the exchange method to exhibit a correlation signal with a low rate of 
overproduction of equal-sized fragment events. In fact, in the BOB calcula- 
tion equal-sized fragment events are produced early in the collision when the 
system goes through the spinodal region of the phase diagram of the nuclear 
matter. Then, subsequent steps like coalescence to form heavier fragments, 
thermal desexcitation and detection filtering, coupled with the finite size of 
the system, deeply blur the original picture of the partitions [24]. Thus, it 
seems difficult to use BOB calculation as a reference to assess the sensitivity 
of the various correlation methods. 

We note that, applied to the same Xe-|-Sn data analysed in this work, the IPM 
method shows an enhancement of equal-sized fragment production at 32, 39 
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and 45 MeV/n [16], in contradiction with our results. 



6 Conclusions 

We have developed a two-step analysis to address the issue of charge correla- 
tions in nuclear collisions. First, we have shown, by using a simulation which 
randomly generates fragments, how the extrapolation to low values of the 
AZ distribution as performed in [14] may induce a spurious enhancement of 
the correlation function at small AZ. We have underlined the difficulties to 
make a quantitative analysis of the charge correlations when the total frag- 
ment charge (Zhound) is not conserved in building the uncorrelated yield Y'. 
Then, to avoid such a problem, we have discussed a new method to build the 
yield Y' by mixing events under the constraints of the conservation of the 
total charge of the fragments and the invariance of charge distributions for 
a given fragment multiplicity. This method has been successfully tested with 
a random break-up simulation and is sensitive to any statistically significant 
deviation from a pure uncorrelated fragmentation process. Finally, it has been 
applied in multifragmentation samples of the ^^^Xe-|-""*Sn central collisions at 
energies between 32 and 50 MeV/n. All charge correlation functions at low 
AZ are equal to unity within two standard deviations and we conclude that 
no significant enhancement of equal-sized fragments has been observed. These 
results are in agreement with those obtained in the reactions Xe-|-Cu and 
Ar-|-Au at 50 MeV/n [13] but at variance with other analyses of the Xe-f-Sn 
data using different estimations of the uncorrelated yield Y' [14,16]. In view of 
these contradictory results, caution should be taken before giving any definite 
conclusion on a signature of spinodal decomposition in nuclear systems by 
charge correlation measurements. 

The experiments ^^^Xe-|-"°*Sn were performed at GANIL by the INDRA col- 
laboration. 
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Table 1 

Number of events in the bin AZ < 1 for each multiphcity and their sum for the 
random break-up calculation and the Xe+Sn "single source" selection at four ener- 
gies. 





Break-up 


32 MeV/n 


39 MeV/n 


45 MeV/n 


50 MeV/n 


M=3 


506 


20 


22 


45 


85 


M=4 


382 


16 


46 


72 


194 


M=5 


100 


5 


15 


53 


145 


M=6 


25 


1 


11 


39 


83 


Sum 


1013 


42 


94 


209 


507 


NT 


659260 


32963 


28088 


25870 


26033 


Sum/NT 


0.154% 


0.127% 


0.335% 


0.808% 


1.948% 



The two last rows give respectively the total number of events (NT) and the per- 
centage of events in the bin AZ < 1. The break-up calculation with 20 times the 
statistics has to be compared with the Xe-|-Sn data at 32 MeV/n. 
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Fig. 1. The first column shows the Zbound distributions for Xe+Sn central collision 
data at 32 MeV/n for each fragment muliplicity M = 3 to 6 (rows). The following 
two columns superimpose Z and AZ distributions of the experimental data (black 
symbols) and the random break-up simulation (lines). The last column is similar 
to the third one with a various range: < AZ < 6. Vertical error bars on the 
data points are statistical. The calculated distributions, performed with a statistics 
20 times larger than the data, have been normalised. Error bars for the simulated 
distributions are small and have not been plotted. 
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Fig. 2. Ratio cl) clS 3j function of AZ between the random break-up simulation yield 
and the calculation of all fragment partitions, both constrained by the 32 MeV/n 
Xe+Sn experimental Z^ound distributions, for each multiplicity M = 3 to 6 and 
their sum (rows); b) similar to a) but as a function of < Z > for the first bin in 
AZ. Error bars are statistical. 
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Fig. 3. Charge correlation analysis with random break-up calculations: a) evolution 
of R{AZ, < Z >) as a function of AZ for four (M,< Z >) pairs; b) evolution of 
R{AZ, Zhound) as a function of AZ for six Zf,ound values, each one containing six 
atomic numbers (scale factors are indicated in parentheses to avoid superimposi- 
tion). The straight lines result from exponential fits in the range 2 < AZ < 8. 
Vertical error bars are statistical. 
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Fig. 4. Fragment charge correlations R at AZ < 1 as a function of Zf,ound with 
random break-up calculations and "background" evaluations of [14] (broken line). 
The calculations are derived from the Xe+Sn Abound distributions at 32, 39, 45 and 
50 MeV/n. Vertical error bars are statistical. Horizontal bars indicate the width of 
the Zi^ound bin. 
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Fig. 5. Evolution of the charge correlations R = Y/Y' by using the exchange mixing 
method as a function of AZ with the random break-up calculations for a) each 
multiplicity from M = 3 to M = 6 and Af = 3, 4, 5, 6; b) same as a) with twice the 
number of events in the bin AZ < 1. Error bars are statistical. 
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Y= random break-up Y= random break-up Y= Xe+Sn 32 MeV/n 

Y'= exchange method Y'= exact I' (|^„„j 1"= exchange method 
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Fig. 6. Deviations from 1 of charge correlation functions R{AZ, < Z >) = Y/Y', 
with multiphcities M = 3, 4, 5 and 6, divided by the statistical error (numerator 
only). The numerator Y corresponds a,b) to the random break-up simulation; c) to 
the 32 MeV/n Xe+Sn data. The denominator Y' is determined by a,c) the exchange 
method; b) the exact calculation of partitions. Normal distribution (0,1) has been 
superimposed on every histogram. 



20 



32MeV/u 39MeV/u 45 MeV/u 50 MeV/u 



1.5 
1 
0.5 

1.5 
1 
0.5 

1.5 

1 

II 

^ 0.5 

1.5 
1 
0.5 

1.5 



M=3 



I .... I .... I 



.••'••••ffj 



>4 



I I I I I 1 1 
M=4 



I .... I .... I 



M=5 







M=6 



AW 



0.5 



M=3,'4,5,6- 



-*~*«*xrT -o, 



I .... I .... I 







- U-n^li 



.lit 



5 10 15 5 10 15 5 10 15 5 10 15 

AZ 

Fig. 7. Experimental charge correlations R versus AZ for Xe+Sn central collisions, 
from left to right as a function of bombarding energy (32, 39, 45 and 50 MeV/n) 
and from top to bottom as a function of the fragment multiplicity, M = 3 to 6. The 
bottom row corresponds to the sum over all multiplicities. Dotted lines indicate 
R = I. Error bars are statistical. 
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Fig. 8. Deviations from of the quantity {R—l)/a as a function of AZ for simulated 
break-up events (first column) and experimental Xe+Sn data from 32 to 50 MeV/n 
bombarding energies (columns 2 to 5), for multiplicities M = 3 to 6 (rows 1 to 4). 
The last row corresponds to the sum over all multiplicities. The dotted lines indicate 
la and 2a deviations. 
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Fig. 9. Experimental charge correlations R in the bin IS.Z < 1 for Xe+Sn central 
collisions, from left to right as a function of bombarding energy (32, 39, 45 and 50 
MeV/n) and from top to bottom as a function of the fragment multiplicity, M = 3 
to 6. The bottom row corresponds to the sum over all multiplicities. All data points 
with more than one event have been plotted. Dotted lines indicate R = 1. Error 
bars are statistical. 



23 



